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Description of Research and Results: 

Analysis of existing infrared spectra from JPL's Mark I High Speed 
Interferometer ( HSI ) balloon flights (1976, 1977) and development of analysis 
techniques applicable to similar data from HSI Mark II balloon flights (1979) 
and from the ATMOS experiment (''Dacelab III). Specific techniques under 
development include line-by-line simulation of the spectra to aid in the 
identification of absorbing gases, simultaneous retrieval of pressure and 
temperature profiles using carefully chosen pairs of CO^ absorption lines, 
and the use of these pressures and temperatures in the retrieval of gas con- 
centration profiles for many absorbing species. A search for new absorption 
features was also carried out, and special attention was given to identification 
of absorbing gases in spectral bandpass regions to be measured by the Halogen 
Occultation Experiment (HALOE). 



Future Plans: 


The temperature and pressure inversion and pas concentration inversion 
algorithms will be further developed and refined, and will be used to retrieve 
profiles from the Ma^k I data and from the Mark II data to become available 
in late 1979. The search for new absorption signatures in the Mark I and Mark 
II spectra vi.M continue, leading to the determination of actual amounts or 
upper limits f >r the stratospheric concentrations of CH-jCJ., Hj,CO and 

OCS. More efficient techniques for line-by-line simulation of tae spectra will 


also be explored. 
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I. INTRODUCTION 


The primary purpose of the proposed research has been to analyze 
existing high-resolution infrared spectra measured by JPL's Mark I High Speed 
Interferometer (HSI) in balloon-borne solar occultation experiments (1976 and 
1977)* and to develop analysis techniques applicable to similar data from the 
Mark II instrument to be flown on a balloon in mid-1979 and data from ATMOS, 
which is now scheduled to fly on Spacelab III (1982). The specific objectives 
of the current study Include identification of major absorbers in the Mark I 
spectra, simulation of the absorption spectra using a line-by-line method for 
comparison to the measured spectra, development of inversion techniques to 
obtain pressure and temperature profiles from CO^ absorption lines, and the 
use of these pressures and temperatures in the simultaneous retrieval of various 
gas concentrations, focusing primarily on the gases and spectral regions impor- 
tant for HALOE (Halogen Occultation Experiment). Unknown features in the existing 
spectra were also to be identified and studied in more detail, with specific 
attention being given to the spectral region of the band of CH^Cfc. In the 
following pages of this report, the above-mentioned objectives will be dis- 
cussed individually, describing both the current status of the research and 
the work planned for the future in each of these areas. 
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II. INVFSTIGATIONr. 


1 . Identification of “inlor Absorbers 

Selected wavelength regions of HSI Mark I spectra from two different 
balloon flights (Texas, 1976 nnd Australia, 1977) are available in graphical form, and 
the complete spectra are stored in digital form on magnetic tape. With the help of 
the Atlas of Infrared Absorption Lines (NASA CR-2925) published at NASA Langley Re- 
search Center, as well as additional data on methane (CK^) provided by Dr. Toth at 
JPL, features belonging to the major absorbing gn3es (H o 0, C0„, 0,, CH^, N ? 0) 
are readily identified in the graphical spectra. A spectral plotting program, 
derived from one developed by R. Norton and B. Stedry at JPL, allows the 
display of spectral data with a scale and resolution chosen by the user 
(Fig. l). Coupled with the line-parameter plotting program used to generate 
the Atlas (Fig. 2), the spectral plotting program facilitates the identifica- 
tion of a large number of absorption features in the Mark I spectra. 

2. Simulation of the Spectra. 

A multi-gas version of the line-by-line model developed by Drayson 
(1966) has been used to simulate several regions of the HSI Mark I spectra 
(Fig. 3). The model has been modified to include the effects of the real 
instrumental line shape for the interferometer, and a fnster routine to 
calculate the Voigt line shape (Drayson, 1976) has been Included. Far.ameters 
in the model can be easily changed in order to simulate the high-resolution 
Mark II spectra when they become available. The agreement between the simu- 
lated spectra find the measured spectra is quite good. The remaining discrep- 
ancies are almost entirely due to the amounts of gases assumed in the model. 
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filter effects present in the measured spectra, and, in some spectral regions, 
missing or uncertain line parameters input into the model. The computational 
efficiency and speed of the model, running on the ('DC 6600 computer, have been 
improved, but the total cost of the simultations remains high. An alternative 
and more efficient simulation technique, based on the Four ier-Trans form method 
of Dr. Mankin at NCAR, which runs on the STAR computer at LoRC, has been 
successfully used in LIMS data analysis at this center. This FFT model will 
be tested and evaluated for use in simulating HSI and ATMOS spectra. 

3. Pressure and Temperature Sensing Using COg Line s 

In order to obtain gas concentration profiles in a limb-viewing 
experiment, it is desirable to get temperature and pressure profiles directly 
from spectral data rather than relying on climatological data and tangent 
heights calculated from instrument pointing information. We have begun develop- 
ing an iterative technique for simultaneous retrieval of pressure and tempera- 
ture profiles using a pair of C0^ absorption lines having different lower- 
state energies, but appearing in the same spectrum. Investigation of the 
temperature and pressure dependence of C0 o line strengths has indicated that 
lines having lower-state energies on the order of 0 -800 cm”^ are best suited 
for pressure-sensing, while those with lower-state energies greater than 
1000 cm * are most strongly temperature-dependent. A technique for temperature 
sensing using isolated CO,, lines in the HSI spectra has already been developed 
by Toth (1977), and pressures may be retrieved using a curve-of-growth method 
(Park et^ al_ 1979) f provided that equivalent widths of the individual C0 o lines 
or a group of lines can be determined with sufficient accuracy. A numerical 
algorithm where the retrieved temperature and pressure profiles are used to 
correct each other in successive iterations has recently been written and has 
been ter'ed using simulated equivalent width values. 
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Figures U and 5 show the expected error in temperature and pressure for 
a single atmospheric layer, retrieved using the above-mentioned algorithm. As 
predicted from analysis of line-strength temperature and pressure dependence, the 
highest accuracies in retrieved temperatures are achieved for CO,, lines having lower- 
state energies greater than 1000 cm For the case shown in Fig. 5 where we assume 
that the temperature profile is known perfectly (AT ■ 0), the error in retrieved 
pressure doe3 not appear to vary with the lower-state energy of the CO^ line chosen 
for pressure sensing. When uncertainty AT is included in the temperature profile, 
the pressure uncertainty Ap remains small for C0^ lines having lower-state energies 
less than 600 cm ^ . Accuracies for both pressure and temperature retrievals 
deteriorate when the absorption line is very weak (absorption less than lOt) or 
nearly saturated (absorption greater than 95?). The error analyses in Figs. 1* and 
5 were carried out for a rather low atmospheric pressure, 0.7 mb, which corresponds 
to an altitude of 50 km. We have found that the accuracy of both the pressure and 
temperature retrievals for higher atmospheric pressures is considerably better than 
the results shown in the figures. 

Testing and further development of the pressure and temperature retrieval 
algorithm will continue, making use of equivalent widths of C0 o lines measured in 
the HSI Mark I spectra. There is also under development a program which will auto- 
matically scan the HSI spectra stored on magnetic tape, find the centers of spectral 
lines and calculate their equivalent widths. The line-finding and equivalent width 
routines are based on the DFRPTH program written by H. Delouis at the Centre National 
de la Recherche Scientifique, Orsay. Eventually these routines will be incorporated, 
along with the atlas and spectral ratio and plotting programs, into a software 
library for analysis of HSI and ATMOS data at LaRC. 

. Simultaneous Retrieval of Gas Concentrations 

Several methods for retrieval of gas concentrations hav been studied, 
among them the direct curve-of-growth method, which requires accurate uetermi- 
nation of equivalent widths, iterative simulation and fitting of line profiles. 
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and the equivalent-width-ratio method used by Fanner and Paper (1977). An 
algorithm for inversion of a gas concentration profile using the single-line 
curve of growth method combined with Drayson's "onion peeling" technique has 
been written and is currently being tested. Further development and testing 
cf gas-concentration retrieval methods will proceed in parallel with the work 
on the temperature and pressure retrievals and automatic calculation of equi- 
valent widths. In addition, selection of absorption lines that can be used 
for most accurate retrieval of concentrations of stratospheric minor gases will 
be carried out. 

5. Support for HALOE 

Work done as part of the current investigation has contributed to 
the support of HALOE, aiding in identification of the best spectral regions 
in which to detect each of the gases of interest (HC£, HF, CH^, NO, h'^O, and 
COg). Detailed analysis and line identification in each of the chosen spectral 
regions remains to be done, particularly in regard to the identification of 
solar lines. Line-by-line simulations have already been carried out for each 
of these regions, and the pressure, temperature and gas concentration retrieval 
algorithms under development in this investigation will be applied first to 
the HALOE spectral regions. The gas concentration profiles determined from 
the HSI spectra will be used to refine previous sensitivity analyses for HALOE. 
Retrieval of HC£ mixing-ratio profiles will receive high priority because the 
HALOE brassboard instrument for the HC£ channel will be test-flown on the CV990 
in the summer of 1979. Since both the ATMOS and HALOE instruments are now 
scheduled to fly on Spacelab III, there will be a unique opportunity for com- 
parison of gas concentration profiles determined from the two experiments. 
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6. Line Identification and rearch for New Signatures 

An extensive analysis and search for new signatures has been carried 
out for the spectral region between 3055 and 3075 cm” 1 . Attention 

was focused on this region to determine if absorption features belonging to 
the fundamental of methyl chloride (CH^CA) were detectable in the HSI Mark I 
spectra. By comparison of the real spectra with Atlas line data and simulated 
spectra, there was found a broad, weak absorption feature extending between 
3062 and 3072 cijT\ which appears in the Texas spectra but not in the Australian 
data. Several species including H^O, C0 2 , 0^ and CH^ have been ruled out as 
the cause of this feature. Simulations of CH^Cl absorption in this region, 
assuming its stratospheric mixing ratios to be one to ten times the values 
calculated by Crutzen et al (1978), seem to indicate that methyl chloride 
alone could not account for the broad envelope found in the Texas spectra. 

The 3Mm absorption spectra of other halogenated hydrocarbons, including methyl 
chloroform (CH^CCfc^), are now being studied and compared to the Mark I data 
in an effort to make a definite identification of the broad feature. Upper 
limits on the stratospheric concentrations of CH^Cfc, CH^CCi,^ and other halo- 
genated hydrocarbons will be calculated after the identification of the 
broad feature is made. 

Preliminary surveys of the Mark I spectra to search for the signa- 
tures of Freon-lL (CF^), formaldehyde (HgCO) and carbonyl sulfide (OCS) have 
also been carried out. Since the strong and absorption bands of CF^ lie 
outside the bandpass of the filters used in the Mark I balloon flights, we 
cannot determine anything about CF^ concentrations from this data. Formaldehyde 
and carbonyl sulfide have strong absorption bands at 3.5Um and 5vm (see Fig. 1), 
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respectively, which have not been detected in the Mark I spectra through 
visual inspection. However, careful analysis of the spectra in these regions 
and simulations of the spectra based on accurate laboratory data may reveal 
weak absorption by OCS or H^CO in these two regions and allow the determina- 
tion of gas concentrations or upper limits for these stratospheric trace 
species. 
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III. rUMMARY 

Considerable p^gress has been made this year in a number of the 
research areas specified in the objectives of this study. Identification of 
absorption features belonging to major infrared-active atmospheric constituents 
in any region of the ’-lark I spectra is accomplished readily by the use of both 
the Atlas and spectral plotting programs. The line-by-line simulation model operat- 
ing on the CDC 6600 computer has been made as efficient as possible, und the less 
costly FFT model running on the STAR computer may be used for simulation of 
HSI data in the near future. Algorithms for the combined retrieval of pressure 
and temperature profiles using CO^ lines and for the retrieval of gas concen- 
trations have been outlined, and early versions of numerical routines for these 
retrievals have been tested. A new absorption feature in the 3.3pm region 
has been found in the Texas data, but is not yet identified. Preliminary 
surveys of the spectra have also been carried out in preparation for the 
determination of upper limits on the stratospheric concentrations of CH^Ci., 

CH^CCfc^, H^CO and OCS. The work on pressure, temperature and gas concentra- 
tion retrievals will be emphasized in the coming year, in support of HALOE 


flight experiments. 
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Figure 2. 
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Figure 5. 


LIST 07 FIGURES 

Plot of a portion of three spectra from the May 1976 HSI Mark I 
balloon flight from Palestine, Texas. The runs 200U, 2000 and 
198k refer to mean tangent heights of 23 km, 33 km, and in excess 
of Uo km, respectively. 

Atlas plot of line data used in the simulation shown in Fig. 3. 

Simulation of Run 200b for the same spectral interval as in the 
previous figure. 

Expected error AT for temperature retrieved for a single layer 
using one CO^ absorption line, assuming that pressure profile is 
known exactly (Ap = 0). E" is the lower-state energy of the line 
in cm~* . 

Expected error Ap for pressure retrieved for a single layer 
using one COg absorption line, assuming that the temperature 
profile is known exactly (AT 0). 
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